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We compare the 7-ray photon flux variability of northern blazars in the Fermi /LAT First Source 
Catalog with 37 GHz radio flux density curves from the Metsahovi quasar monitoring program. 
We find that the relationship between simultaneous millimeter (mm) flux density and 7-ray photon 
flux is different for different types of blazars. The flux relation between the two bands is positively 
correlated for quasars and does no exist for BLLacs. Furthermore, we find that the levels of 7-ray 
emission in high states depend on the phase of the high frequency radio flare, with the brightest 
7-ray events coinciding with the initial stages of a mm flare. The mean observed delay from the 
beginning of a mm flare to the peak of the 7-ray emission is about 70 days, which places the average 
location of the 7-ray production at or downstream of the radio core. We discuss alternative scenarios 
for the production of 7-rays at distances of parsecs along the length of the jet . 



I. INTRODUCTION 

Several studies, based on the first year of 
Fermi/LAT operations, have shown that: (i) the 7- 
ray and the averaged radio flux densities are signifi- 
cantly correlated [e.g. 1 and (ii) blazars detected at 
7-rays are more likely to have larger Doppler factors 
and larger apparent opening angles than those not 
detected by LAT [e.g. [2]. This observational evidence 
strongly suggests that radio and 7-ray emission have 
a co-spatial origin. To locate and identify the region 
where the bulk of 7-ray emission is produced, and 
to provide details about its connection to the radio 
jet, an analysis of simultaneous radio and 7-ray light 
curves is essential. 

However, we highlight two caveats about the inter- 
pretation of radio/gamma correlation analyses, which 
have often been interpreted too simplistically. First, it 
is well-known that there is usually a considerable de- 
lay between mm and cm radio flares. Thus, although 
cm- flares would tend to peak after the 7-ray flares, 
mm-flares would show shorter delays or possibly even 
peak before the 7-rays. The very important second 
caveat is that a correlation analysis tends to measure 
the distance between the peaks, especially if the flares 
have different timescales (as the radio and the 7-ray 
flares tend to have ). However, a radio flare starts to 
grow a considerable time before it peaks. The begin- 
ning of a millimeter flare coincides with the ejection 
of a new VLBI component from the radio core [3]. 
This is the epoch that must be compared with the 
7-ray flaring, not the epoch of the radio flare maxi- 
mum. The crucial question is whether a 7-ray flare 



occurs before the beginning of a mm- flare, or after it; 
in the former case the 7-rays originate upstream of 
the radio core (the beginning of the radio jet), in the 
latter, they originate downstream of the radio core, 
presumably from the same disturbances that produce 
the radio outbursts. 

In this proceeding we summarize our recent re- 
sults [4], obtained after combining the finely sam- 
pled 37 GHz Metsahovi light curves and the monthly 
binned 7-ray light curves provided by the Fermi/LAT 
First Source Catalog [1FGL, [5]. By using a radio 
flare decomposition method, we estimate the begin- 
ning epochs of millimeter flares and their phases dur- 
ing 7-ray flaring events to establish the true temporal 
sequence between 7-ray and radio flaring. 



II. THE 37 GHZ AND 7-RAY LIGHT CURVES 

The Metsahovi quasar monitoring program cur- 
rently includes about 250 AGN at 37 GHz. From 
them, we selected a sample of sources that fulfill the 
following criteria: (i) well-sampled light curves during 
the period 2007-2010, covering the 1FGL period, (ii) 
a firm association with the 1FGL catalogue, and (iii) 
a 7-ray monthly light curve during the 1FGL period 
that is significantly different from a flat one. 

Our final sample consists of 60 sources classified ac- 
cording to their optical spectral type as highly po- 
larized quasars (HPQ, 22), low polarization quasars 
(LPQ, 5), quasars without any optical polarization 
data (QSO, 15), BL Lac type objects (BLO, 17), and 
radio galaxies (GAL, 1). The sample of sources used 
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FIG. 1: Monthly flux-flux relation for the combined sample of 60 northern AGN. The different types of sources are 
symbol coded as shown in the legend, the correlations coefficients are shown only when the significance is > 99.9%. 



in this work is listed in Table 1 of [4 



III. THE S 3 7GHz - S 7 CORRELATION 

To compare radio flux densities and 7-ray photon 
fluxes, monthly binned radio light curves were created 
from the Metsahovi light curves. The time bins are 
the same as in the 1FGL flux history curves, allowing 
us to compare simultaneously 7-ray photon flux and 
radio flux-density variations with the time resolution 
of a month. 

Figure 1 shows that simultaneous measurements 
at the two bands appear to be positively corre- 
lated. However, we find significant differences between 
quasars and BL Lacs, which we describe below. By 
applying the Spearman's rank correlation test, two 
very clear results emerge. First, there is a significant 
positive correlation between the 7-ray photon flux and 
the 37 GHz flux density for quasars, while the BL Lac 
fluxes are not correlated. Second, the strength and 



the significance of the correlations is different for each 
type of quasar. 

The photon flux - flux density correlation for 
quasars is absent for QSOs, significant for LPQs, and 
very significant for HPQs. Such a dependence on the 
degree of optical polarization may arise naturally if 
the polarization indicates the viewing angle of the 
jet, with sources with high optical polarization hav- 
ing their jets oriented closest to our line of sight (see 
Figure 7 in [6]). The dependence of the flux - flux 
relation on optical polarization agrees with previous 
results, where it has been shown that the brightest 
7-ray emitters have preferentially smaller viewing an- 
gles and consequently higher Doppler boosting factors 
[e.g., [2]. Since 7-ray fluxes and radio flux densities 
are significantly correlated for sources where the rel- 
ativists jet is aligned close to our line of sight, this 
implies that there is a strong coupling between the 
radio and the 7-ray emission mechanisms. That the 
correlation is seen on monthly timescales further indi- 
cates a cospatial origin in quasar- type blazars. 
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FIG. 2: The recent flux history at 37 GHz (filled circles) and 7-rays (triangles) of 3C 345. The total flux density has 
been decomposed using individual exponential flares. The filled triangle represents the most significant peak in the 7-ray 
flux density during the 1FGL period. The dashed vertical line is drawn to highlight the relation between the peak in the 
7-ray light curve and the ongoing millimeter flare. 



IV. CONNECTIONS BETWEEN ONGOING 
FLARES AND HIGH STATES OF 
GAMMA-RAY EMISSION 



We have decomposed the mm light curves into indi- 
vidual exponential flares (see Figure 2), each of which 
corresponds to a new disturbance created in the jet 
and is often detectable as a new VLBI component [3] . 
We further calculate the phase of the mm-flare when 
the most prominent maxima in the 1FGL light curves 
occurred. Our analysis shows that the most signifi- 
cant 7-ray flux peaks tend to occur when a mm-flare 
is either rising or peaking. This indicates that the 
strong 7-ray flares are produced in the same distur- 
bances that produce the mm flares. 

Using our data, we can estimate the time delay be- 
tween the time when the mm flare starts and when the 
7-rays peak for each source. We define the beginning 



of a mm flare to be: 

irnrn j_mm /i \ 

L — L max T 5 \ L ) 

where t 7 ^^. is the time of the mm flare peak and r 
is the variability timescale. In other words, we define 
the beginning of the flare as the epoch when its flux 
is 1/e of the maximum flux, S^™ 171 ) = Sm e ax . For 
each source we estimate the time delay between the 
time of mm-flare onset (green circle in Figure 2) and 
when the 7-ray peak occurs (filled triangle in Figure 
2). The observed time delay has a distribution cen- 
tered around 70 days with the onset of the mm-flare 
preceding the 7-ray peak. 

After converting the time delays to linear distances 
from the region where the mm-outburst begins (i.e. 
the radio-core) to the region of the 7-ray production, 
our estimates lead us to conclude that in our sample 
the average location of the 7-ray emission region is 
about 7 parsecs downstream the radio-core. 
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FIG. 3: Sketch of the inner parts of an AGN. Recent 
results [131 114j have shown that a significant amount of 
broad-line emission is driven by optical continuum radia- 
tion produced downstream of the radio core, this strongly 
suggests the presence of an additional component of the 
BLR located at parsec scales. In this scenario, an outflow- 
ing BLR may serve as a source of seed-photons to produce 
7-rays via EC models far away from the central engine [3]. 



V. SUMMARY AND DISCUSSION 

The results presented in [4] strongly indicate that at 
least for the strongest 7-rays the production sites are 
downstream or within the radio core [7], well outside 
the BLR at distances of several parsecs or even tens of 
parsecs from the black hole and the accretion disk. A 
number of papers based on Fermi data have reached 
similar conclusions [e.g. [8]. 

In the current AGN paradigm, it is widely believed 
that 7-ray are produced via Inverse Compton (IC) 
mechanisms in the relativistic jet, and the seed pho- 
tons for the IC process may be provided by the jet syn- 
chrotron emission (Synchrotron Self Compton, SSC) 
as well as by external photon sources (External Comp- 
ton, EC) such as an accretion disk, the broad-line re- 



gion (BLR) or the hot dusty torus. At the distance 
of 7 parsecs from the radio core - keeping in mind 
that the radio core itself is often at a considerable dis- 
tance from the black hole [e.g. [9] - the only sources of 
seed photons for the IC processes are the jet itself and 
the dusty torus. However, there is growing evidence 
[13j[14] that the BLR might extend to much larger dis- 
tances than given by virial estimates (Rblr < 1 pc). 
One possibility is that the jet drags a part of the BLR 
whit it, see Figure 3 for a sketch of an outflowing BLR 
among the other inner AGN constituents. 

While single zone SSC has failed to reproduce the 
observed 7-rays [10] and so far there are only a couple 
of blazars with firm detections of a dusty torus [TT| fT2] . 
a tentative idea to test is whether an outflowing BLR 
can serve as a source of external photons to produce 7- 
rays, even at distances of parsecs downstream of the 
radio- core. In this scenario, the strong 7-ray events 
are produced in the same disturbance that produces 
the radio outburst by upscattering external photons 
provided by an outflowing BLR. 

The most effective way to explore any of the above 
scenarios (and others) is by modeling simultaneous, 
well-sampled spectral energy distributions (SEDs) 
[TSJ [16] following a multizone modeling approach [17] . 
Furthermore, recent results have suggested that the 
more massive the black hole is, the faster and the 
more luminous jet it produces [18 . Therefore, a reli- 
able estimate of the black hole masses is an essential 
input to theoretical models of both the shape and the 
variability of blazars SEDs. 
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